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Abstract 

Radiative transfer equation in an accretion disk wind is examined analytically and numerically under 
the plane-parallel approximation in the subrelativistic regime of (v/c) 1 , where v is the wind vertical velocity. 
Emergent intensity is analytically obtained for the case of a large optical depth, where the flow speed and 
the source function arc almost constant. The usual limb-darkening effect, which depends on the direction 
cosine at the zero-optical depth surface, does not appear, since the source function is constant. Because 
of the vertical motion of winds, however, the emergent intensity exhibits the velocity- dependent limb- 
darkening effect, which comes from the Doppler and aberration effects. Radiative moments and emergent 
intensity are also numerically obtained. When the flow speed is small (v < 0.1c), the radiative structure 
resembles to that of the static atmosphere, where the source function is proportional to the optical depth, 
and the usual limb-darkening effect exists. When the flow speed becomes large, on the other hand, the 
flow speed attains the constant terminal one, and the velocity-dependent limb-darkening effect appears. 
We thus carefully treat and estimate the wind luminosity and limb-darkening effect, when wc observe an 
accretion disk wind. 
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stars 



1. Introduction 

Accretion disks are now widely believed to be energy 
sources in various active phenomena in the universe: 
in protoplanctary nebulae around young stellar objects 
(YSOs), in cataclysmic variables (CVs) and supersoft X- 
ray sources (SSXSs), in galactic X-ray binaries and mi- 
croquasars (//QSOs), and in active galaxies (ANGs) and 
quasars (QSOs). Accretion-disk models have been ex- 
tensively studies during these three decades (see Kato 
et al. 1998 for a review). Besides the traditional stan- 
dard disk by Shakura and Sunayev (1973), new type disks, 
such as advection-dominated accretion flows (ADAF) or 
radiatively-inefficient accretion flows (RIAF) for the very 
small mass-accretion rate (e.g., Narayan, Yi 1994), and 
supercritical accretion disks or so-called slim disks for the 
very large mass-accretion rate (e.g., Abramowicz et al. 
1988). 

Accretion disk winds have been also extensively exam- 
ined in relation to astrophysical jets and outflows: in bipo- 
lar outflows from YSOs, in mass outflows from CVs and 
SSXSs, in relativistic jets from /iQSOs, AGNs, QSOs, and 
in gamma-ray bursts (GRBs). In particular, intense radi- 
ation fields of luminous supercritical accretion disks may 
be responsible for relativistic jets from super-Eddington 
sources, such as luminous //QSOs, GRS 1915+105 and 
SS 433, luminous QSOs, 3C 273, and energetic GRBs (see, 
e.g., Fukue 2004 for references). 

In such circumstances, radiative transfer in accretion 
disk winds as well as accretion disks becomes more and 
more important. 



Radiative transfer in the standard disk has been inves- 
tigated in relation to the structure of a static disk atmo- 
sphere and the spectral energy distribution from the disk 
surface (e.g., Meyer, Mcyer-Hofmeister 1982; Cannizzo, 
Wheeler 1984). Furthermore, gray and non-gray mod- 
els of accretion disks were constructed under numerical 
treatments (Kfiz and Hubeny 1986; Shaviv and Wchrsc 
1986; Adam et al. 1988; Mineshige, Wood 1990; Ross et 
al. 1992; Shimura and Takahara 1993; Hubeny, Hubeny 
1997, 1998; Hubeny et al. 2000, 2001; Davis et al. 2005; 
Hui et al. 2005) and under analytical ones (Hubeny 1990; 
Artcmova et al. 1996; Fukue, Akizuki 2006a). 

Radiative transfer in the accretion disk wind, on the 
other hand, has not been well considered both in the non- 
relativistic and relativistic regimes. Recently, radiative 
transfer in a moving disk atmosphere was firstly investi- 
gated in the subrelativistic regime (Fukue 2005a, 2006a), 
and in the relativistic regime (Fukue 2005b, 2006b; Fukue, 
Akizuki 2006b). In contrast to the static atmosphere, in 
the moving atmosphere the boundary condition at the 
surface of zero optical depth should be modified (Fukue 
2005a, b). Moreover, the usual Eddington approximation 
violates in the highly relativistic flow (Fukue 2005b; see 
also Turolla, Nobili 1988; Nobili et al. 1991; Turolla et 
al. 1995; Dullemond 1999), and the velocity-dependent 
variable Eddington factor was proposed (Fukue 2006b for 
a plane-parallel case; Akizuki, Fukue 2007 for a spherical 
case). 

Radiation hydrodynamical (RHD) simulations were 
also performed for radiation-dominated supercritical disks 
with winds by several researchers (Eggum et al. 1985, 
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1988; Okudaet al. 1997, 2005; Okuda, Fujita 2000; Okuda 
2002; Ohsuga et al. 2005; Ohsuga 2006). In these current 
studies of RHD simulations for disks and winds, they were 
done in the subrclativistic regime up to the order of (w/c) 1 , 
using the moment formalism and the flux-limited diffusion 
(FLD) approximation (Levermore, Pomraning 1981). The 
flux-limited diffusion method provides good approxima- 
tions to the exact solutions but only if they are derived 
from transfer equations in which terms of the order of 
(v/c) 2 or higher have been retained (Yin, Miller 1995). 

Radiative transfer problems on accretion disks and 
winds are not well understood yet, in particular for the 
relativistic cases. Hence, in order clarify the physics, in 
addition to RHD simulations, we must treat the simplified 
problem in the analytical way. 

In this paper, we thus examine radiative transfer in 
the accretion disk wind, which is assumed to blow off 
from the luminous disk in the vertical direction (plane- 
parallel approximation), and analytically and numerically 
obtain the flow solutions for the case without internal 
heating. In the previous studies, radiative flows in the 
vertically moving atmosphere were solved in the subrela- 
tivistic regime (Fukue 2005a, 2006a), and in the relativis- 
tic regime (Fukuc 2005b, 2006b; Fukue, Akizuki 2006b), 
where only radiative moments were obtained. In the 
present paper, we further obtain the radiation intensity 
as well as radiative moments. 

In the next section we describe the basic equations. In 
section 3, we show analytical solutions, while wc present 
numerical solutions in section 4. The final section is de- 
voted to concluding remarks. 

2. Basic Equations 

Let us suppose a luminous flat disk, deep inside which 
gravitational or nuclear energy is released via viscous heat- 
ing or other processes. The radiation energy is trans- 
ported in the vertical direction, and the disk gas, itself, 
also moves in the vertical direction as a disk wind due 
to the action of radiation pressure (i.e., plane-parallel ap- 
proximation). For simplicity, in the present paper, the 
radiation field is considered to be sufficiently intense that 
both the gravitational field of, e.g., the central object and 
the gas pressure can be ignored. We also assume the gray 
approximation, where the opacities do not depend on the 
frequency. As for the order of the flow velocity v, we con- 
sider the subrclativistic regime, where the terms of the 
first order of (v/c) are retained, in order to take account 
of radiation drag. 

The radiative transfer equations are given in several 
literatures (Chandrasekhar 1960; Mihalas 1970; Rybicki, 
Lightman 1979; Mihalas, Mihalas 1984; Shu 1991; Kato 
et al. 1998). For the plane-parallel geometry in the ver- 
tical direction (z), the radiation hydrodynamic equations 
are described as follows (Kato et al. 1998). It should be 
noted that the basic equations below are the same as those 
given in Fukuc (2005a), except for the transfer equation 
for the radiation intensity. 

For matter, the continuity equation is 



pv = J (= const.), (1) 

where p is the gas density, v the vertical velocity, and J 
the mass-loss rate per unit area. The equation of motion 
is 



dv_ _ Kabs + 

dz c 



[F-{E + P)v]. 



(2) 



where K a bs and K sca are the absorption and scattering 
opacities (gray), which are defined in the comoving (fluid) 
frame, E the radiation energy density, F the radiative 
flux, and P the radiation pressure in the vertical direc- 
tion, which are measured in the fixed (laboratory) frame. 
In a gas-prcssurcless approximation, the energy equation 
is reduced to 



= q-p[j- CKabs-E + Kabs 



2Fv 



(3) 



where q + is the heating and j is the emissivity, which is 
measured in the comoving frame. 

For radiation fields, the frequency-integrated transfer 
equation, the zeroth moment equation, and the first mo- 
ment equation become, respectively, 



dl 

P~r=P 
dz 



(l + 3/? M ) 



An 



(K abs + K SC a)(l-/?/i)/ 



+ ■ 



-in 



{(l + 3/3p,)cE-2F/3} 



dF 

dz 
dP 

dz 



p (j - cK ahs E + n &hs Ff3 - K sca Fp) , 



P 



[j(3 - (K abs + K sca )F 

-K 3hs cP/3 + K Bca (cE + cP)j3]., 



(4) 
(5) 

(0) 



where /i is cos 6*, 6 being the polar angle, and (3 = v/c. We 
further adopt the Eddington approximation in the comov- 
ing frame, which is transformed into 



1 



4F 



P=-E+--[3 (7) 

3 6 c 

in the fixed frame (Kato et al. 1998). Here, the transfer 
equation (4) is corrected to the order of {v/c) 1 (Kato et 
al. 1998). 

Eliminating j with the help of equation (3), and intro- 
ducing the optical depth by 



dr- 



■(k, 



abs 



i)pdz, 



(8) 



we can rearrange the basic equations up to the order of 
{v/c) 1 as 



„2rW 



c z J^- = -(F-AcP(3), 
dr 

dl 1 a+ 



(9) 



"3F = 



47T («abs + «sca)p 



(1+3/?//) 

1 



+ (l-0li)I-(l + 30ii) — (cE-2F0), (10) 

47T 



dF 
dP 



( K abs + K sca,)P 



■F0, 



( K abs + K sca,)P 



/3 + F-4cPf3, 



(11) 
(12) 
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J 



,dz 



1 



-c/3. 



(13) 



Finally, integrating the sum of equations (9) and (12) 
gives the momentum flux conservation in the present ap- 
proximation. 

c 2 J(3 + cP = cP - j- — /3dr = cP , (14) 

J ( K abs + K sca )p 

when there is no heating (q + = 0). In addition, the sub- 
script means the value at some reference position (i.e., 
the wind base). Similarly from equations (9) and (11) we 
have the energy flux conservation, 



\jv 2 + F = F - 

.1 (, K abs T K sc& )p 



dr = F , 



(15) 



when there is no heating (q + = 0). Here, the first term 
on the left-hand side is eventually dropped, although we 
retain it here to clarify the physical meanings. 

We solve equations (9), (14), (15), and (10) for appro- 
priate boundary conditions, and we obtaine analytic solu- 
tions. 

As for the boundary conditions at the wind base of r = 
To and at the wind top of r = 0, we impose the following 
conditions. 

At the wind base on the disk surface with an arbitrary 
optical depth To, the flow velocity j3 is zero, the radiative 
flux is Fq (which is a measure of the strength of radiation 
field), and the radiation pressure is Pq (which connects 
with the radiation pressure gradient), where the subscript 
denotes the values at the wind base. 

At the wind top, on the other hand, as already pointed 
out in Fukue (2005b), the usual boundary conditions for 
the static atmosphere cannot be used for the present ra- 
diative wind, which moves with velocity at the order of the 
speed of light. Namely, the radiation field just above the 
wind top changes when the gas itself does move upward, 
since the direciton and intensity of radiation change due 
to relativistic aberration and Doppler effect (cf. Kato et 
al. 1998; Fukue 2000). If a flat infinite plane with surface 
intensity 7 S in the comoving frame is not static, but mov- 
ing upward with a speed v s (= c/3 s , and the corresponding 
Lorentz factor is 7 S ), where the subscript s denotes the 
values at the surface, then, just above the surface, the 
radiation energy density E a , the radiative flux F s , and 
the radiation pressure P s measured in the inertial frame 
become, respectively, 



cE s = 27t7, 



F s = 2tt7 



cP s = 2tt7 



3 7 2 + 3 7s ii s 



3 7s 2 + 8 7s u s + 3m 2 
6 

7 2 + 3 7s w s + 3m 2 



(16) 
(17) 
(18) 



where u s (~j s v s /c) is the flow four velocity at the surface 
(Fukue 2005b). As a result, we have the boundary condi- 
tion at the wind top within the present approximation, 



cP 
F m 



2 + 6/3 s + 6/? s 2 

3 + 8/3 s + 3/3 2 



-A- 



(19) 



This is the consistent boundary condition under the 
present subrclativistic regime up to the order of (y/c) 1 , 
although Fukue (2005a) have approximately used the 
boundary condition, cP s /F s = 2/3, for a static atmosphere. 
In general, at the wind top of t = 0, the boundary condi- 
tion (19) is not satisfied. Hence, for given parameters, we 
adjust and obtain the mass-loss rate J as an eigen value, 
so as to satisfy the boundary condition (19). 

As already stated, for the internal heating in the wind, 
we assume that there is no heating source (q + = 0), al- 
though it is straightfoward to extend the model to the 
case with heating. On the other hand, at the wind base 
on the luminous disk, there is assumed to be a uniform 
source of 7o . 

3. Analytical Solutions 

In Fukue (2005a), analytical solutions for radiative mo- 
ments were already derived. For the completeness, in 
this section, we first recalculate them under the consistent 
boundary condition. Using the analytical expressions for 
moments, we then calculate the radiative intensity, which 
we wish to know in the present paper. 

3.1. Flow Velocity and Radiative Moments 

Under the present subrclativistic regime, equation (15) 
means that the radiative flux F is conserved: 

F = F = F S . (20) 

Using the boundary conditions at the wind base (v = 0), 
equation (14) is expressed as 

Jv + P^Pq. (21) 

Hence, equation (9) becomes 



"i~« 



47» , 
that can be analytically solved to yield 



4P 



1 



4P 



(T - T ) 



(22) 



(23) 



Thus, the radiative flow from the luminous disk without 
heating is expressed in terms of the boundary values and 
the mass-loss rate. In addition, the flow velocity v s at the 
wind top (t = 0) is 

4Pn 




1-e cJ 



T() 



(24) 



Using the boundary condition at the wind top, we 
further impose a condition on the values at boundaries. 
Inserting boundary values (19) into momentum equation 
(21), using equation (23), we have the following relation: 

/ 47}, 



P 



cJ 

47b 



Fs 



18cPn 



1-e cJ 



TO 



(25) 



That is, for given To and Po at the wind base, the mass- 
loss rate J is determined in units of F s /c 2 , as an eigen 
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value. Compared with Fukue (2005a), the second term in 
the first parentheses on the right-hand side is an additional 
one, which appears due to the present corrected boundary 
condition (19). 

As already stated in Fukuc (2005a), the mass-loss rate 
increases as the initial radiation pressure increases, while 
the flow terminal speed increases as the initial radiation 
pressure and the loaded mass decrease. 

Moreover, as easily shown from equation (25), in order 
for the flow to exist, the radiation pressure P$ at the flow 
base is restricted in some range, 



2 + V6 cP 2 



(26) 



which is slightly modified from that given in Fukue 
(2005a), due to the corrected boundary condition. At the 
upper limit of cPo/F s = 2/3 + r , the loaded mass diverges 
and the flow terminal speed becomes zero. On the other 
hand, at the lower limit of cP$/F s = (2 + \/6)/6, the pres- 
sure gradient vanishes, the loaded mass becomes zero and 
the terminal speed approaches the saturation speed of 



/3oo = 



-6 + 3^6 



0.337, 



(27) 



where the radiative flux F is balanced by the radiation 
drag APqv with the boundary condition (19). 

In figure 1 we show analytical solutions, the flow veloc- 
ity v (solid curves) in units of c, the radiation pressure 
P (dashed ones) in units of F s /c, and the source function 
S (chain-dotted ones) in units of F s as a function of the 
optical depth r for several values of Pq at the wind base 
in a few cases of tq. 

When the initial radiation pressure Pq at the wind base 
is large, the pressure gradient between the wind base and 
the wind top is also large. As a result, the loaded mass J 
also becomes large, but the flow final speed v s is small due 
to momentum conservation (14). When the initial radia- 
tion pressure Pq is small, on the other hand, the pressure 
gradient becomes small, and the loaded mass is also small, 
but the flow final speed becomes large. In the latter case, 
the source function 5" becomes almost constant. 

Within the present approximation of (v/c) 1 with the 
boundary condition (19), the flow final speed saturates at 
0.337 c. 

3.2. Emergent Intensity 

Now, we turn to the transfer equation (10), or 

dl , „ . , „ , 1 



M^ = (l-/W-(l + 3/3/i)^(c£- 



2F0), 



(28) 



when there is no heating. 

Under the present approximation up to the order of 
(v/c) 1 , F is constant, but E (P) and (3 are generally func- 
tions of r. If, however, the radiation field is sufficiently 
intense, the wind flow quickly saturates; the radiative 
flux, F, is balanced by the radiation drag force, (E + P)v, 
and the wind speed reaches the saturation terminal one, 
F/(E + P). In such a case, the wind speed and the radi- 
ation quantities are almost constant. In other words, the 
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Fig. 1. Flow velocity v (solid curves) in units of c, radiation 
pressure P (dashed ones) in units of F B /c, and source function 
S in units of F s (chain-dotted ones) as a function of the optical 
depth t for several values of Pq at the wind base in a few cases 
of to: (a) to = 1 and (b) to = 10. From top to bottom of v 
and from bottom to top of P and S, the values of Po are 0.75, 
1, 1.2, 1.4 in (a), and 0.75, 1, 2, 5 in (b). 

source function, the second term on the right-hand side of 
equation (28), is almost constant (see figure 1), since in 
the present case the source function S is expressed as 

s= c_E-2F l= 3cP-6F l 

47T 47T 

Here, we thus assume that the wind velocity is constant, 
and analytically integrate the transfer equation (28). 

When the wind velocity becomes a constant terminal 
one, using equations (23), (21), and (7), we have 



0=A 



4cP n 



cP = cP - c 2 J(5 =-F s + —^- 



cE = 'icP - 4F(3 ■■ 



Hence, 



3 
2F S 

AFl 2 



1 F 2 

__ S 

18 cP ' 
5F[_ 
6cP Q - 



16 



cE - 2F/3 = 2F s --^ r = 2F s - —F s f3. 
S cPq 3 



(30) 

(31) 
(32) 

(33) 
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That is, the source function measured in the fixed frame 
slightly decreases due to the effect of the relativistic mo- 
tion. 

Under the above situations, we can now integrate the 
radiative transfer equation (28), similar to Fukue and 
Akizuki (2006). After several partial integrations, we ob- 
tain both an outward intensity I(r,fj.) (/i > 0) and an in- 
ward intensity I(t, — /x) as 

cE - 2FI3 1 + 3/3/i 



I(t,h) 



J(t, -aO 



4tt 1-Ph 

+/(r ,/i)e i ^ ii(r " ro) 
cE - 2F/3 1 + 3/3/i 
4^ 1-Afy 



1 



-{T-T ) 



1 



(34) 
(35) 



where I(tq,h) is the boundary value at the wind base on 
the luminous disk. 

In general case with finite optical depth To and uniform 
incident intensity Iq from the disk, the boundary value 
J(to,/x) of the outward intensity I consists of two parts: 

I(to,h)=I + I{t ,-h), (36) 

where Iq (= F s /tt) is the uniform incident intensity and 
I(tq,— fx) is the inward intensity from the backside of the 
disk beyond the midplanc. Determining I(tq,~ fi) from 
equation (35), we finally obtain the outward intensity as 

cE - 2F[3 1 + 3/3/x 



Hr,fi) = 



1 



4tt 1-PfJ, 



3F S (2 16 8 



S(t-2t ) 



1-e" 



1-/3M 



(t-2t ) 



4 i^i 
3 



(r-r„) 



(37) 



where we have used F s = ttIq . 

Finally, the emergent intensity 1(0, [i) emitted from the 
wind top becomes 



J ^) = S 



r 



2 16 „ 8 \ / 1-/J*. 



2r 



3F s /2 16 8 ' 



for large r . (38) 



Under the present approximation, where the source 
function is constant, the usual limb darkening does not 
appear: e.g., (3F s /4tt)(2/3 + t) for the Milne-Eddington 
solution. Due to the Doppler and aberration effects orig- 
inating from the vertical motion of winds, however, the 
emergent intensity (38) depends on the wind velocity as 
well as the direction cosine. This is the velocity-dependent 
limb-darkening effect. 

In figure 2, the emergent intensity 1(0, fj,) normalized by 
the isotropic value I (= F s /tt) is shown for several values of 
(3 as a function of fj,. Although the present approximation 
may be valid for j3 < 0.1, we show the cases for /3 < 0.3 in 
order to stress the velocity-dependency. 




Fig. 2. Normalized emergent intensity as a function of /j, 
for the case without heating. The numbers attached on each 
curve are values of f} of wind velocity. The dashed straight line 
is for the usual Milne- Eddinton solution for the plane-parallel 
case. 

As is easily seen in figure 2, as the velocity becomes 
large, the limb-darkening effect becomes prominant. That 
is, the emergent intensity increases in the poleward direc- 
tion, while it decreases in the edgeward direction. As a 
result, a wind luminosity would be overestimated by a 
pole-on observer and underestimated by an edge-on ob- 
server, when we observe an optically-thick accretion disk 
wind. 

It is interesting that the intensity does not change in 
the direction of fi = 2/3. 

4. Numerical Solutions 

In this section, we numerically solve equations (9) and 
(12) [or (14)] with constant F, and further solve equation 
(10) [or (28)] for several cases, and compare the results 
with that of analytical solutions to check the accuracy 
and limitations of analytical solutions. 

4--1- Flow Velocity and Radiative Moments 

In figure 3 we show several numerical solutions, the flow 
velocity v (thick solid curves) in units of c and the radia- 
tion pressure P (thick dashed curves) in units of F s /c as 
a function of the optical depth t for several values of Pq 
at the wind base in a few cases of tq . Corresponding an- 
alytical solutions are shown by thin curves, although the 
mass-loss rates are slightly different so as to satisfy the 
upper boundary condition (19). 

As is seen in figure 3, in both analytical and numerical 
cases, when the initial pressure Pq and the loaded mass 
are large, the flow speed is small (v < 0.1c). In such a case, 
the radiative structure resembles to that of the static at- 
mosphere, where the source function is proportional to the 
optical depth. When the initial pressure and the loaded 
mass decrease, on the other hand, the flow speed becomes 
large to attain the saturation one. In such a case, the 
source function is almost constant. 

In addition, for the same parameters the analytical 
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Fig. 3. Flow velocity v (thick solid curves) in units of c and 
radiation pressure P (thick dashed curves) in units of F B /c 
as a function of the optical depth r for several values of Po 
at the wind base in a few cases of to: (a) to = 1 and (b) 
to = 10. From top to bottom of v and from bottom to top of 
P, the values of Po are 0.75 and 1.4 in (a), and 1 and 5 in (b). 
Corresponding analytical solutions are shown by thin curves. 

and numerical solutions arc slightly different, although wc 
have solved the same basic equations. This is understood 
as follows. Our basic equations are up to the order of 
(v/c) , and therefore have the accuracy of the same or- 
der. In deriving the analytical solutions, we have dropped 
the terms of order of (v/c) 2 , while we have retained those 
terms in solving the numerical solutions. Hence, the ana- 
lytical and numerical solutions should be same within the 
accuracy of (v/c) . 

4-2. Emergent Intensity 

After obtaining the numerical solutions for v and P (i.e., 
E), we can further numerically integrate the transfer equa- 
tion (28) for many angles /i as a function of r under ap- 
propriate boundary conditions and given initial conditions 
of tq and Po- As boundary conditions, wc set a uniform 
source of Iq at the wind base, while there is no incident 
intensity at the wind top: / + (to,/i) = Iq + /~(to,— /i) and 
J - (to,— /i) =0. We use meshes of 100 for /i and 500 for r. 
Finally, the emergent intensity 7(0, /i) emitted from the 
wind top is numerically obtained as a function of angle 



Fig. 4. Normalized emergent intensity as a function of fi 
for the case without heating. The numbers attached on each 
curve are values of /3 S of wind terminal velocity. The ini- 
tial optical depth tq is 1 and the initial pressure Po is 1.4 
(/3s = 0.150) and 0.75 (/3s = 0.337). 

for given tq and Po- Examples of the results are shown in 
figure 4. 

In figure 4, the emergent intensity obtained numerically 
is shown for the initial optical depth of tq = 1, and for 
several values of (3 Sl that is attached on each curve. When 
the terminal speed is small, as already stated, the wind 
structure resembles to that of the static atmosphere. In 
the present case, however, the optical depth of wind is 
finite. As a result, at around /x~ 1 the emergent intensity 
is almost constant with the values of Jo (i.e., the peaking 
effect diminishes), while the usual limb-darkening effect 
recovers for small /x, where the line-of-sight length is long. 

When the terminal speed is large, on the other hand, the 
velocity-dependent limb-darkening effect appears again 
for the numerical solutions. That is to say, as stated in 
figure 2, due to the Dopplcr and aberration effects orig- 
inating from the vertical motion of winds, the emergent 
intensity depends on the wind velocity as well as the di- 
rection cosine. In particular, the emergent intensity is 
enhanced and increases as \i increases. 

By the way, in figure 4 the emergent intensity entirely 
decreases for all /u, as the velocity increases. Apparently, 
this seems to be curious, because the relativistic boosts 
work in the regime with large velocity. The reason is also 
the relativistic effect. The present subrelativistic flow up 
to the order of (v/c) , the momentum and energy of ra- 
diation fields accelerate the flow. Although the flux F is 
conserved within the present order of (v/c) 1 , the radia- 
tion intensity diminishes due to the interaction between 
the field and the flow [the 2P/3 term in the second term 
on the right-hand side of equation (28)]. In other words, 
the source term effectively decreases. As a result, the 
emergent intensity entirely decreases, compared with the 
non-relativistic case. 

We here briefly check the consistency of the numerical 
results. In order to obtain the radiation intensity / from 
equation (28), we use the radiation energy density £ mom 
and the radiative flux F mom (— F s ) obtained from the mo- 
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Fig. 5. Ratios between the quantities obtained from the mo- 
ment equations and those obtained from the intensity. Solid 
curves and dashed ones denote the energy density and the 
flux, respectively. In the small velocity case of /3 S = 0.150 
(thick curves), both quantities are almost consistent. In the 
large velocity case of /3 B = 0.337 (thin curves), on the other 
hand, the consistency becomes worse. In addition, chain- 
dotted curves mean the three times Eddington factor for the 
quantities obtained from the intensity. 

ment equations (14). Once the radiation intensity I(r,[i) 
is obtained, we can calculate the radiation energy den- 
sity -Esim and the radiative flux F s i m by the definition. If 
the quantities obtained from the moment equations and 
those obtained from the intensity coincide each other, the 
solutions are consistent, and vice versa. 

In figure 5 we plot the ratios B s im/fimom (solid curves) 
and F s i m / F mom (dashed ones) for the small velocity case 
of (3 S = 0.150 (thick curves) and for the large velocity case 
of P s = 0.337 (thin curves). In addition, we also plot 
the quantities (.E s i m + 4-F S i m v)/(3P s i m ), the three times 
Eddington factor, by chain-dotted curves. 

As is seen from figure 5, in the case of small veloc- 
ity, where the structure resembles to that of the static 
case, the ratios are almost unity and the solutions are 
consistent. In the small velocity case, furthermore, the 
Eddington factor is almost 1/3. In the case of large veloc- 
ity, on the other hand, the ratios differ from unity. The 
main reason of this discrepancy is the relativistic effect. 
As already stated, the radiation intensity diminishes due 
to the interaction between the field and the flow. As a 
result, the intensity I(j, fj,) as well as the emergent inten- 
sity 1(0,^) entirely decreases. Thus, F s [ m becomes smaller 
than F mom . Although the reason that the energy density 
ratio increases is less clear, it seems to be the relativis- 
tic peaking effect (Doppler boost and relativistic abbera- 
tion). Indeed, the Eddington factor in the large velocity 
case is somewhat larger than unity. This would be also 
the relativistic peaking effect. Anyway, as the velocity 
becomes large, the consistency becomes worse, and the 
present treatment under the order of (v/c) 1 would be in- 
valid. 



In this paper we have examined the radiative trans- 
fer problem in an accretion disk wind under the plane- 
parallel approximation in the subrelativistic regime of 
(v/c) 1 . The flow velocity, the radiation pressure distribu- 
tion, and other quantities are analytically and numerically 
solved as a function of the optical depth for the case with- 
out heating. Furthermore, the emergent intensity from 
the wind top with zero optical depth is also analytically 
obtained under the assumption of constant flow velocity, 
while numerically calculated. 

When the source function is constant, te usual limb- 
darkening effect does not appear. However, the emergent 
intensity exhibits the velocity-dependent limb-darkening 
effect, which originates from the Doppler and aberration 
effects, associating with the vertical subrelativistic motion 
of winds. As a result, a wind luminosity would be over- 
estimated by a pole-on observer and underestimated by 
an edge-on observer, when we observe an optically thick 
accretion disk wind. 

It should be noted that the apparent optical depth in the 
rclativistically moving media. Abramowicz et al. (1991) 
pointed out that the optical depth in the relativistic flow 
decreases as 7(1 — /3/j)t toward the downstream direction, 
due to the Doppler and aberration effects. Inspecting 
equation (28) or solution (38), we find that, in the present 
subrelativistic flow, the optical depth r is apparently re- 
placed by (1 — 0/J,)t. This is just consistent with the results 
by Abramowicz et al. (1991) within the order of (v/c) 1 . 

In the present modelized plane-parallel geometry, the 
gas density at the wind top is finite, since the wind veloc- 
ity there is finite. Rigorously speaking, this situation is 
inconsistent with the definition of the wind top, where the 
optical depth is effectively zero. In more realistic cases, 
the disk wind would geometrically diverge, and the den- 
sity quickly drops. As a result, the surface of a moving 
photosphere exists, even if the wind terminal velocity is fi- 
nite. The present model mimics such a realistic disk wind, 
where the velocity-dependent limb darkening effect would 
also appear. In other words, as long as the acceleration 
of wind takes place in a small height, and then the wind 
diverges, the present assumption would be valid. 

The radiative transfer problem investigated in the 
present paper must be quite fundamental problems for ac- 
cretion disk physics and astrophysical jet formation. In 
this paper, for simplicity, we only considered the subrela- 
tivistic case of (v/c) 1 . The fully relativistic case of (v/c) 2 
will be explored in the future. 
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